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Abstract 

The light cone QCD sum rules are employed to calculate the 



7riViV(1535) coupling g n NN*- We study the two point correlation 
function of two nucleon currents sandwiched between the vacuum 
and the pion state. The contribution from the excited states and the 



^ . continuum is subtracted cleanly through the double Borel transform 

with respect to the nucleon and iV(1535) momenta, pf, p\ = (p—q) 2 - 



Our calculation shows that the 7riViV(1535) coupling is strongly 



! suppressed. 

PACS Indices: 13.75.Gx; 14.20.Gk; 14.40.Aq; 13.75. Cs; 12.38.Lg 
Quantum Chromodynamics (QCD) is asymptotically free and its high energy behav- 
ior has been tested to one-loop accuracy. On the other hand, the low-energy behavior 
has become a very active research field in the past years. Various hadronic resonances 
act as suitable labs for exploring the nonperturbative QCD. The inner structure of nu- 
cleon and mesons and their interactions is of central importance in nuclear and particle 
physics. Internationally there are a number of experimental collaborations, like TJNAL 
(former CEBAF), COSY, ELS A (Bonn), MAMI (Mainz) and Spring8 (Japan), which will 
extensively study the excitation of higher nucleon resonances. 

Among various baryon resonances, negative parity resonance iV*(1535) is particularly 
interesting, which dominates the rj meson photo- or electro-production on a nucleon. 
The branching ratio for the decay N* — > r]N is comparable with that for N* — > ttN. 
Considering the phase space difference and using the experimental decay width of N* ||, 
we get QrjNN* = 2 and g n NN* = 0.7. The latter value is in strong contrast with the pion 
nucleon coupling g n NN = 13.4. Thus arises the question: why is the coupling g n NN* so 
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small compared with g n NN^ 

Whether the coupling g n NN* is strongly suppressed is under heated debate in literature. 
In a recent extended coupled channel analysis of tcN scattering, the Jiilich group used 
9rjNN* = 1-94 and g^NN* < 0.12 j4|. Jido, Oka and Hosaka suggested in a recent letter 
that irNN* coupling is strongly suppressed as a consequence of chiral symmetry 0] . Their 
argument was based on a pion-nucleon correlator of two baryon interpolating fields. The 
chiral transformation properties of their interpolating fields then imply that the correlator 
is purely proportional to the tensor structure 75, with no piece of the formes, which is the 
relevant structure for ttNN* coupling. Based on the above observation they claimed that 
the coupling ttNN* vanishes. The above argument was criticized by Birse |J. He pointed 
out that the P75 piece of the correlator is a sum of all possible pion-baryon couplings 
that can contribute. Hence the absence of a ^75 piece is a statement about the particular 
combination of the pion-baryon coupling and the subtraction terms that correspond to 
the chosen interpolating fields. It does not imply that the physical ttNN* coupling is 
suppressed. With an interpolating field with covariant derivative for the iV* Kim and Lee 
used QCD sum rules to estimate g^NN* ~ 1-5 |7|]. But in their analysis the continuum 
and excited states contribution is poorly subtracted and the numerical results depend 
strongly on the value of the quark gluon mix condensate (qg s cr ■ Gq), which renders their 
conclusion less convincing. 

Although it is widely accepted that QCD is the underlying theory of the strong in- 
teraction, the self-interaction of the gluons causes the infrared behavior and the vacuum 
of QCD highly nontrivial. In the typical hadronic scale QCD is nonperturbative which 
makes the first principle calculation unrealistic except the lattice QCD approach, which is 
very computer time consuming. So a quantitative calculation of the 7riV\/V(i535) coupling 
with a tractable and reliable theoretical approach proves valuable. 

Among the various nonperturbative methods, QCD sum rules (QSR) is very useful in 
extracting low-lying hadron masses and coupling constants [|TJ. The light cone QCD sum 
rules (LCQSR) is quite different from the conventional QSR, which is based on the short- 
distance operator product expansion (OPE). The LCQSR is based on the OPE on the 
light cone, which is the expansion over the twists of the operators. The main contribution 
comes from the lowest twist operator. Matrix elements of nonlocal operators sandwiched 
between a hadronic state and the vacuum defines the hadron wave functions. When the 
LCQSR is used to calculate the coupling constant, the double Borel transformation is 
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always invoked so that the excited states and the continuum contribution can be treated 
quite nicely. Moreover, the final sum rule depends only on the value of the wave function 
at a specific point like ip w (uo = |), which is much better known than the whole wave 
function ||. 

LCQSR has been used to derive the couplings of pions with heavy mesons in full QCD 
HI, in the limit of tuq — > oo || and I/ttlq correction [[TTJ, the couplings of pions with 
heavy baryons [II], and various semileptonic decays of heavy mesons |T2[ etc. 



We shall employ the LCQSR to calculate nNN* coupling. The continuum and excited 
states contribution is subtraced cleanly within our approach. 
We start with the two point function 

n(pi,p 2 ,g) = / d A xe^{0\Tr] p {x;s)if n {0;t)\7r + {q)) (1) 

with pi = p, p 2 = p — q and the general nucleon interpolating field without derivatives 
which couples to both positive and negative parity nucleon resonances 

rjp(x; s) = e abc { \u a (x)Cd b (x)} 7 5 w c (x) + s \u a (x)C^d b (x)} u c (x)} (2) 

where a, b, c is the color indices, C = i^lo is the charge conjugation matrix, s is the 
mixing parameter and rj p = ^ p ^7o- For the neutron interpolating field, u d. The Ioffe's 
current 7] p (x; s — — 1) couples strongly to the positive parity nucleon [Q, O], while it was 
found that the current rj p (x; s = 0.8) is optimized for negative parity nucleons and couples 
strongly to JV(1535) 0. 

U(pi,P2,q) has the general form 

n (pi , p 2 , q) = F (pj , pi, q 2 ) g7 5 + Fx (jp\ , pj, q 2 ) 75 + F 2 (p\ , pi, q 2 )p^ 5 + F 3 (pj , p\ , q 2 ) o^p^q v 

(3) 

It was well known that the sum rules derived from the chiral odd tensor structure yield 
better results than those from the chiral even ones in the QSR analysis of the nucleon 
mass and magnetic moment 0, 15]. Most of the QSR analysis of the pion nucleon coupling 
constant deals with the tensor structure 575. It is important to note that the diagonal 
transitions like iV — > N, N* — > N* does not contribute to the tensor structure P75. In 
other words, the function F 2 involves solely with the process iV* — > N and corresponding 
continuum contribution. Based on the above observation we shall focus on the chiral odd 
tensor structure pj 5 throughout. 
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The 7riV7V(1535) coupling constant g n NN* is denned as: 

£"kNN* = 9-kNN*Nt ■ ttN + h.c. (4) 
At the phenomenological level the eq.((l|) can be expressed as: 

n / \ / . x r X N (s)X N *(t) \ N *(s)\ N (t) 

U( Pl ,p 2 ,q) = -(m N +m N *) 9 « NN *{ ¥i _ _ m^ - J^ _ Mlt){pl _ Mf) >^+- • 

(5) 

where we write the structure P75 explicitly only and the continuum contribution is denoted 
by the ellipse. \n(s) is the overlapping amplitude of the interpolating current t]n{x) with 
the nucleon state 

(0|t7(0; s)\N(p)) = \ N (s)u N (p) (6) 

By the operator expansion on the light-cone the matrix element of the nonlocal op- 
erators between the vacuum and pion state defines the two particle pion wave function 
(PWF). Up to twist four the Dirac components of this wave function can be written as 



< n(q)\d(x) W u(0)\V > = t du e luqx (^(u) + x 2 9l (u) + 0(x 4 )) 

<J 

+ ^) f'due^g^u), (7) 



qx Jo 



< 7r(g)|rf(x)i7 5 M(0)|0 > = J7r n / due luqx (p P (u) , (8) 

m u + m d Jo 

— f IT? /"I 

< n(q)\d(x)cr^ 5 u(0)\0 > = i(q^x u - q u x^) n n du e luqx ip a {u) . (9) 

6(m„ + m d ) Jo 



< Tr(q)\d(x)a a 0j 5 g s G IJiU (ux)u(O)\O >= 

if3K[{QnQa9up ~ QvQagtf) - {%%9v a - qugpg^a)} / Van <f 37r (a i )e tqx(oil+va3) , (10) 



< ir(q)\d(x)7^ 5 g s G a(3 (vx)u(0)\0 > 



UUpigav — ) - qaigpti — )1 / Va^^a^e 

LV q ■ x ' v q ■ x n J 

+U^(q a x l 3-q l3 x a ) [ Va mi {a i )e iqx ^+ va ^ 
q ■ x J 



D iqx{a\+votz) 



and 



< ■K{q)\d{x)^ IJi g s G ol p{vx)u{Q)\Q > 



4 



g • x ' v g ■ x ' 

+if«-^-(q a Xfi - g^a) / Pa i ^||(a i )e^ ai+Ma ) . (12) 
g • x j 

The operator G a p is the dual of G a p: G a/ 3 = \e a i35pG sp ; Vcti is defined as Voii = 
doiida.2doiz8(l — a.\ — a<i — 03). Due to the choice of the gauge x^A^ix) = 0, the path- 
ordered gauge factor Pexp (ig s Jq 1 dux^ A^iux)) has been omitted. The coefficient in front 

of the r.h.s. of eqs. (|8|), (|9]) can be written in terms of the light quark condensate < uu> 

m 2 2 

using the PCAC relation: fi^ = - — = — — < uu >. 

m u + m d f* 

The PWF (p % (u) is associated with the leading twist two operator, gi(u) and ^(w) 
correspond to twist four operators, and (fp(u) and (p a (u) to twist three ones. The function 



(ps-x is of twist three, while all the PWFs appearing in eqs.([TT|), (|12l) are of twist four. The 
PWFs (p(xi,/j,) (n is the renormalization point) describe the distribution in longitudinal 
momenta inside the pion, the parameters x, (J2i x % = 1) representing the fractions of the 
longitudinal momentum carried by the quark, the antiquark and gluon. 

The wave function normalizations immediately follow from the definitions (|7|)-([T2|): 
f^du <f n (u) = Jodu ip a {u) = 1, fodu gi{u) = S 2 /12, JVa^^ai) = JVaiip\\(ai) = 
0, / T>aiip±(ai) = — j T>Oii{p\\(oii) = 5 2 /3, with the parameter 5 defined by the matrix 
element: < Ti{q)\dg s G a ^ a u\Q >= i5 2 

The full light quark propagator with both perturbative term and contribution from 
vacuum fields reads 

iS(x) = {0\T[q(x),q(0m 
. x (qq) x 2 



2vr 2 x 4 12 192 

ig. 



(qg s a ■ Gq) (13) 



/ du{^a ■ Giux) - Uu^G^(ux) lv \ + 

JO X 2 X 1 



167T 2 

where we have introduced x = x^ . 

We neglect the four particle component of the pion wave function and express (|l|) with 
the PWFs. After Fourier transformation and making double Borel transformation with 
the variables p\ and p\ the single-pole terms in (||) are eliminated and finally we arrive 
at: 



(m N + m N *)[\ N (s)\ N *(t) - AAr*(s)A A r(t)]5f 7r ATAr*e 



r " N 1 N* ■> 
77T) 



h (1 +*)(! + tWMM^f^) + -£-(1 + ,)(! + t)^iK)M 4 /i(^) 



192tt 2V Jy '™ uy J ^M 2 ' 16tt 2V yv uy J1 "M 2 



, )0 _ 2 [7(1 + st) + 3(a + t)]<feMM 4 /i(-^) + (1 + a + t - 3st)a^' a (u )M 2 / (^) 



fn 



4608tt 2 



[l + s)(l + t)(g" s G 2 )^(u )M 2 f [ 



so ■ 
M 2 ' 



+ 



1927T 5 



t)amlfj, n ipp(u ) - 



6912vr 2 



(5 + 3s + 3£ — llstyamlfinip'^Uo) 



2304tt 2 



[19(1 + st) + 7 ( s + t)}(g 2 s G 2 }g 2 (u ) 



[l + s)(l + t)(g 2 G 2 )g[(u ) (1,4) 



11527T 2 

where ^ = 1.65GeV, U = 132MeV, (qq) = -(225MeV) 3 , (qg s a ■ Gq) = m 2 (qq), m 2 = 

n , 

0.8GeV 2 , a = — (2-a) 2 (qq). f n (x) = 1 — e~ x tt is the factor used to subtract the 

k=0 

continuum, which is modeled by the dispersion integral in the region s\, s 2 > sq, so is the 



continuum threshold. u 



M l n M 2 



M\Ml 
M'(+M% ' 



M 2 , Mf are the Borel parameters, and 



¥4M 



du 



etc. 



In obtaining ([14]) we have used integration by parts to absorb the factor (q ■ x), which 
leads to the derivatives of PWFs in the above formula. Moreover we have used the double 

m _ [M 2 f- n e-^b{u - u ). 



Borel transformation formula: £>i ^ £>■ A ~ 



2 p| [m 2 — (1— u)p\— up 2 ]" 

We note that the twist four PWFs ip±(ea), ip\\(aii), ( f±(°^i) an d do not contribute 

to the chiral odd tensor structures, which was first observed in [1161. Moreover the twist 



three PWF tp^ appears in the combination f^p^-K{qq)q 2 ■ In the physical limit q 2 



mz 



0, its contribution is negligible, which is in contrast with the sum rule for g^NN il6fl ) where 
the twist three PWF (p^ plays an important role. 

It is a common practice to adopt the Ioffe current for the ground-state nucleon [0], i.e., 



letting s — — 1. Now the equation (14") has a simple form: 

(m N + m N *)[\ N (s)\ N *(t) - \ N *(s)\ N (t)}girNN*e 



-) 



£)<? 2 (u Q )M 4 /i 



M 2 ' 



7^2*°/^^ M M 2 f ( ) 



3456tt 2 



1 + 7t)aml^ip' a (uo) 



1927T 2 

/tt 



1927T 



(1 + t)amlfi n ip P (u ) 
Al-W s G 2 )g 2 (u ) 



(15) 



where £ = 0.8. 

The various parameters which we adopt are a 



0.546 GeV 3 , (g 2 G 2 ) = 0.474 GeV 4 



rrir 



0.8 GeV 2 , s = 3.2GeV 2 , X N {s = -1) = 0.013GeV 3 , \ N *{s = 0.8) = 0.27GeV 3 with 



the formulas in JTJ. Moreover Ajv(s = -1) > Ajv(£ = 0.8), \ N *(t = 0.8) > \ N *(s = -1), 
so it is reasonable to discard the term Ajv(£ = 0.8)Aat*(s = —1) in the eq. flT5|). We use 
the model PWFs presented in [§] to make the numerical analysis. 

Our sum rule ([TBI) is asymmetric with the Borel parameters M 2 and Mf due to the 
significant mass difference of iV and iV(1535). It is natural to let Mf = 2m 2 sf /3, Mf = 



G 



2m 2 N *f3 1 where (3 is a scale factor ranging from 1.0 to 2.0. In this case we have M 2 = 
1.28/?GeV 2 , M = 0.27, g 2 (u ) = -0.03GeV 2 , ¥p(u ) = 0.66 and <^(w ) = 2.63 at the 
scale fj, = lGeV. 

The dependence of g^NN* on the continuum threshold and the Borel parameter is weak. 
It is very stable with reasonable variations of so and M 2 as can be seen in FIG 1. Finally 
we have 



We have included the uncertainty due to the variation of the continuum threshold and 
the Borel parameter j3 in ([OJ). In other words, only the errors arising from numerical 
analysis of the sum rule (p~5|) are considered. Other sources of uncertainty include: (1) 
the truncation of OPE on the light cone and keeping only the lowest twist operators. 
For example the four particle component of PWF is discarded explicitly; (2) the inherent 
uncertainty due to the detailed shape of PWFs in different models; (3) throwing away 
the term A;v(£ = 0.8)Aat» (s = —1) in the eq. (|T5|); (4) the continuum model etc. In the 
present case the major error comes from the uncertainty of PWFs since our final sum 
rule (|I~5D depends both on the value of PWFs and their derivatives at uq = 0.27. Luckily 
around Uq = 0.27 different model PWFs have a shape roughly consistent with each other. 
So a more conservative estimate is to enlarge the error by a factor of two. Now we have 



We want to point out that one should not be too serious about the specific number. 
What's important is the fact that \g n NN*\ "C IflW/wl- Our result is consistent with the 
conclusions in Refs. [§], [5]. 

In summary we have used LCQSR to calculate the irNN* coupling. The continuum 
and the excited states contribution is subtracted rather cleanly through the double Borel 
transformation with respect to the nucleon and iV(1535) momenta, p 2 , p 2 = {p — q) 2 - Our 
result shows explicitly the suppression of g^NN*- 
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9nNN* = (-)(0.08±0.03) . 



(16) 



9nNN* = (-)(0.08±0.06) . 



(17) 
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Figure Captions 

FIG 1. The sum rule for g^NN* as a function of the scale parameter (3 with the continuum 
threshold s = 3.4, 3.2, 3.0 GeV 2 using the PWFs in [§]. 
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